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Abstract 

We study the cosmic no-hair in the presence of spin-2 matter, i.e. in bimetric gravity. We obtain 
stable de Sitter solutions with the cosmological constant in the physical sector and find an evidence 
that the cosmic no-hair is correct. In the presence of the other cosmological constant, there are 
two branches of de Sitter solutions. Under anisotropic perturbations, one of them is always stable 
and there is no violation of the cosmic no-hair at the linear level. The stability of the other 
branch depends on parameters and the cosmic no-hair can be violated in general. Remarkably, the 
bifurcation point of two branches exactly coincides with the Higuchi bound. It turns out that there 
exists a de Sitter solution for which the cosmic no-hair holds at the linear level and the effective 
mass for the anisotropic perturbations is above the Higuchi bound. 
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I. INTRODUCTION 



It is well recognized that the large scale structure of the universe stems from primordial 
fluctuations generated quantum mechanically during inflation. Remarkably, the nature of 
primordial fluctuations is independent of initial conditions. This nice feature can be asso- 
ciated with the conjecture that the initial anisotropy and inhomogeneity rapidly disappear. 
This is called the cosmic no-hair conjecture. The cosmic no-hair is proved in an ideal situa- 
tion [l| . Namely, a homogeneous expanding spacetime with a cosmological constant rapidly 
approaches de Sitter spacetime, i.e., the initial anisotropy decays in a Hubble time, when 
we assume that matter satisfies the strong and dominant energy conditions. In general, 
however, it is not clear whether the cosmic no-hair conjecture is correct or not. In fact, a 
counter example to this conjecture was found |2|. There, spin-1 gauge fields remain during 
inflation and the anisotropy does not necessarily vanish. Moreover, it turned out that anti- 
symmetric tensor fields can also generate the anisotropy 3|. Hence, it is natural to explore 
the possibility that a symmetric spin-2 tensor as matter causes the violation of the cosmic 
no-hair conjecture. Historically, a model of massive spin-2 matter has been proposed as that 
of meson J4], which can be regarded as bimetric gravity consisting of the physical metric 
and the other spin-2 tensor field. In order to treat the spin-2 matter, therefore, we need to 
construct a consistent ghost-free theory of bimetric gravity. Fortunately, this task has been 
accomplished recently p-19]. 

Given a consistent model of spin-2 matter, we can study the cosmic no-hair conjecture. 
There are some reasons that we expect the conjecture can be violated. In the presence of 
spin-2 matter, it is inevitable that gravitons have the mass as the consequence of mixing 
between the physical metric and the other spin-2 tensor field. When we consider massive 
gravitons in an expanding spacetime, the decay time scale of the anisotropy is determined 
by comparing Hubble scale with the effective mass of gravitons. For example, by taking 
the couplings of the physical metric and the spin-2 matter small, the Hubble friction term 
might be dominant compared with the effective mass term in the equation of motion then 
the decay time scale becomes much longer than Hubble time scale. Besides the above one, 
there may be the violation of the energy conditions in the presence of the spin-2 matter 10]. 
Since the energy conditions are assumed in the proof of [l|, it is not apparent whether the 
cosmic no-hair holds or not in bimetric gravity. 
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In this paper, we consider a cosmological constant in bimetric gravity as the limit of 
slow roll inflation. First, we concretely reveal the property of de Sitter solutions in bimetric 
gravity. Then, we investigate the fate of the anisotropy perturbatively. We stress that it is 
important to study the background geometry in detail because the effective mass of gravitons 
can depend on the background geometry. Since, in known cases, the violation of the cosmic 
no-hair appears already at the linear level, we expect the linear analysis reflects the feature 
at the nonlinear level. 

When we consider massive gravitons in de Sitter spacetime, we also need to care about 
the fact that the helicity-0 mode of massive gravitons becomes a ghost when the effective 



mass is below the 
Deser type ghost 



iiguchi bound llMl4j| . Note that this ghost is different from a Boulware- 



15] which is already removed by construction. Since there is no a priori 



reason to forbid the mass of gravitons violating the Higuchi bound, we also check if the 
effective mass satisfies the Higuchi bound. 

We organize the paper as follows. In section Ull we present ghost-free bimetric gravity and 
derive basic equations needed for the analysis. In section IHIl we study the cosmic no-hair in 
bimetric gravity in the presence of the cosmological constant in the physical sector. We find 
that de Sitter solutions are stable and the small anisotropy rapidly decays. In section IIVt we 
introduce the other cosmological constant and investigate the stability of de Sitter solutions 
and the fate of the anisotropy. We also study whether the Higuchi bound is satisfied or not. 
The final section is devoted to the conclusion. In appendix A, we derive a set of equations 
used in the text. 



II. BIMETRIC GRAVITY 



In this section, we introduce bimetric gravity as a model of spin-2 matter and 

provide basic formulae. Historically, after the pioneering work J4], bimetric gravity has been 



studied from time to time 
multi-spin-2 matter 
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17| . The model can be generalized to that of ghost-free 



19]. 



Let us represent the physical metric and the other metric as g^ v and f^ u , respectively. 
Note that we regard the other metric as the spin-2 matter. We consider bimetric gravity 
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with cosmological constants 
M 2 



M„ f . Mi f , . 

S = -f J d 4 xV^9(R[9,u] - 2A 9 ) + -± J d x\f^f{R[f j 

r _ 3 

+m 2 M e 2 / d 4 xv^J> n F n [L£] , 

J n=l 



2A, 



(1) 



where M g and Mf are Planck constants of g^ and f^, and R is the scalar curvature 
constructed from each metric. The interaction terms of the metrics are defined as 



This combination of interaction terms gives no Boulware-Deser ghost |q|. Here, m 2 is a 
coupling constant of the metrics and {a n } n =i t 2,3 are arbitrary constants. We define the 
reduced Planck constant M e as 

1 1 1 

Mf ~ M] + M] ' 



20] when we take the 



where M e is chosen so that m coincides with the Fierz-Pauli mass 
massive gravity limit. Note that we can regard A g as the potential energy of a scalar field 
in the slow roll approximation coupled to the physical metric g^ v as in general relativity. 

In this paper, we consider the simplest case a<i — 1, a,\ = — 0. Then, the action is 
written as 



5 



M 2 r 

-± J d 4 Xy/^(R[g^} - 2A g ) 
+m 2 M 2 J d 4 x^F 2 [L^] , 



Mf 



d^-f\R\f, u ] - 2h f ) 



(2) 



where 

m»] = \{\L} 2 - [L 2 ]) , [L] = , [L 2 ] 

We now present basic equations and derive formulae which will be used in the later 
analysis. 



A. de Sitter solutions in bimetric gravity 



In this subsection, we consider homogeneous and isotropic solutions in bimetric gravity 
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21 



- |25| . We derive equations of motion and show the solutions are de Sitter spacetimes. 



We take the homogeneous and isotropic metric ansatz for g^ u and f^, 

ds 2 = -N 2 (t)dt 2 + e 2ait) [dx 2 + dy 2 + dz 2 } , 



(3) 



and 



ds ' 2 = —M 2 (t)dt 2 + e 2m [dx 2 + dy 2 + dz 2 } , (4) 

respectively. M, N are lapse functions and a, /3 describe the isotropic expansion of each 
metric. Substituting the metric ansatz into the action, we obtain the Lagrangian 

3/3 2 



C 



2„3a 



■lf- NA : 



M 



MA 



+m 2 M 2 Ne 3a [6 - 9e + 3e 2 + 7 (-3 + 3e 



where 
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M 
~N 



(5) 



(6) 



Taking the variation with respect to each variable, we obtain the equations of motion for a 

and (3 



M 



and two constraints 



(£) - M M-m)(-- t ) = o 

' + «l-o 9 )£- 3 (M-Af e )(2- e ) 

GvO = \, + M2-e)(e-l) , 
(§) 2 = A/ + £(l-a 9 )e- 3 (l-e) 

where we normalized parameters and time with M e as follows: 
M 2 m 2 . A 
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A, 



3M e 2 ' 



A 



/ 



3M C 2 ' 



(7) 
(8) 

(9) 
(10) 



-- • (ID 

M e dt 



We notice that a g can take the value in the range < a g < 1 from the definition of M e . The 
detailed derivation can be found in Appendix 

In bimetric gravity, the diagonal part of general coordinate invariance is preserved. Hence, 
the two constraints contain a first class constraint and a second class constraint. Thus, there 
exists a secondary constraint. Now, from ([7j) and (jUJ) (or (jSJ) and ffTUj) ). we can deduce the 
equation 

/3 \ /RipP a'e a \ 

^ =0. (12) 



<(!-« 



Af 



A/ 



The first factor can be taken to be zero. However, this is a special solution and it is known 



that this leads to a pathology 
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-|30|. Hence, we take the following branch 



M = —Ne . (13) 

a 



This is nothing but the condition determining the Lagrange multiplier. From (Q, ffTUj) and 
( TIB"]) , we obtain the secondary constraint 

g{e) — (A/ + £a g )e 3 - 3£a 9 e 2 + [-X g + 2£a g - f (1 - a g )]e + {(1 - a g ) = . (14) 

From the definition of e, e should be positive and hence we should look for the positive roots 
of the algebraic equation g(e) = 0. Since £, a g , X g and Xf are constants, a positive root 
of g(e) = is also a constant which we represent eo- Then, taking the derivative of the 
definition of e, we derive a' = (3' and hence M = Ne^. Now, we take a gauge N — 1 using 
the gauge degree of freedom. Then, we get M = eo = constant. From ((Tj) and (jHJ), we can 
deduce a" = (3" = which can be solved as a = H Q M e t, (3 = H Q M e t + log (eo), where H is 
Hubble scale which is determined from the constraints as 

Hi = A 9 + £a g (2-e )(e -l) 

= A/eg + tfl-a,)^. (15) 

Thus, we obtained two de Sitter spacetimes with the relation = Cq provided that eo 
is a positive root of g(e) = and Hq > holds for eo- 



B. Fate of the anisotropy 

In this subsection, we consider the anisotropy perturbatively and examine how the 
anisotropy evolves. We also derive the effective mass of the massive graviton. 
We take the anisotropic metric ansatz 

ds 2 = -N\t)dt 2 + e 2a{t) [e'^dx 2 + e 2a{t) {dy 2 + dz 2 )] , (16) 

and 

ds' 2 = -M 2 {t)dt 2 + e 2 ^[e^dx 2 + e 2X ®{dy 2 + dz 2 )} , (17) 

where a and A describes the anisotropic expansion of each metric. Here we assume the 
anisotropy is small. Substituting the metric ansatz into the action and dropping the higher 
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order terms, we can derive the quadratic Lagrangian 



5 2 C = M 2 e ia ^- + M|e 3 ^ + m 2 M 2 e Ne 3a [-9e + 3e 2 + 3 7 e]g 2 , (18) 



where we defined the new variable 



M 



q = A — a . 



(19) 



Note that a and A can be regarded as zero modes of gravitons. From the above action, we 
can deduce the equations for a as 



and for A as 



a" + 3H a' - £a 9 e (3 - 2e )q = 



A" + 3H X' + f (1 - a fl )-(3 - 2e )q = 



By taking the difference of (12 ip and (1201 . it is easy to obtain 

1 



g // + 3/f ?' + ^ 



a«eo + (1 - Ofl 



(3 - 2e )q = 



(20) 



(21) 



(22) 



From this equation, we can read off the effective mass of the massive graviton as 

1 



m 



off 



a q e + (1-0,0 



(3 - 2e ) . 



(23) 



Since the effective mass is different from the bare mass £, it is non-trivial if the effective 
mass is less than Hubble scale even if the bare mass is so. By making the combination 
(120]) x l/a 3 + (12TT) x e 2 /(l - a g ), we have 

o A' 



D 3H Q t 



a 



e V e r 



1 — a, 



. 



(24) 



This leads to a conserved quantity 



E = e 3Hot (— + el—— — 

\a g 1 — (lg 



(25) 



which means the mode 



a 



X 1 



1 1 

1 - O n 



(26) 



corresponds to the massless graviton. The existence of the massless mode is a reflection of 
the diagonal general coordinate invariance. From the conservation law (1251) . we see that this 
mode vanishes exponentially fast. 



If we substitute q = e lut into eq. (1221) . we obtain 

q = A exp iu + t + B exp 2o;_t , (27) 

where 

.3H , l~ 9HZ 
u ± = i^± ^m 2 ff - (28) 

and A, B are integral constants. If m 2 ff is negative, q exponentially grows like 



/ /. , , 9# 2 3H \ 
q-Bexp (^|m 2 ff | + _° _ . (29) 

Inversely, if m 2 ff is positive, q exponentially decays. When m 2 ff — > 0, the decay time 
scale r is r = 2/3H . On the other hand, if m 2 ff — < 0, the time scale is evaluated as 



r = (XL 



i -^-Vx + (2 ^ 2 - m - ) - (30) 

Therefore, the decay time scale of the anisotropy r is shorter than Hubble time scale 1/H Q 
for m 2 ff > 2Hfi and the opposite holds for < 2Hq. 



III. DECAY OF THE ANISOTROPY: CASES X f = 

First, we consider the situation Xf = 0. The constant X g can be regarded as the poten- 
tial energy of a scalar field coupled to g^ in the slow roll approximation. We prove that 
there exist a de Sitter solution for X g > and the solution is stable under the anisotropic 
perturbations. We also see that the effective mass of the massive graviton is bounded from 
below m 2 ff > 3Hq. This suggests that the anisotropy rapidly decays in a Hubble time. 

When we take Xf = 0, (TH]) and (IT5|) become 

g(e) = £a g e 3 - 3£a 9 e 2 + [-X g + 2£a g - {(1 - a g )]e + £(1 - a g ) = (31) 

and 

# 2 = A 9 + ^(2-e )(e -l) 

= £(l-a 9 )— • (32) 

From the second line of (132)) . we see that eo should be less than 1 so that H is a real number. 
Then, from the first line of (132|) . A 5 should have a positive lower bound. We assume that X g 
is positive in the following. Then, we obtain 

g(0) = £(1 - a g ) > , #(1) = — A fl < , #(e) ->■ +oo as £-)■ +oo . 
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Thus, there are a positive root smaller than 1 and a positive root larger than 1. The root 
larger than 1 does not satisfy the condition Hq > 0. It turned out that there is a single 
positive root eo in the range < e < 1 where H Q is a real number in the case A 9 is positive. 

Next, we will see that the de Sitter solution derived above is always stable under the 
anisotropic perturbations. Apparently, massless modes rapidly decay in Hubble time scale. 
Then the stability under the perturbations of a and A is determined by the sign of the mass 
term of the perturbation equation ff22l) as mentioned in Sec. Ill Bl or the sign of 

m cff = i \ a 9 e o + (1 - a g)— I ( 3 - 2e o) • 
L e J 

Since the de Sitter solution satisfies < e < 1, m 2 s is positive. Therefore, the de Sitter 
solution is stable under the perturbations of a and A. 

Furthermore, we can prove that m 2 s is bounded from below m 2 s > 3Hq. To show this, 
let us define 



in 



2 \„ A 

eff 



[a fl eg + (l-a fl )](3-2e ) 



»(*>=ir • V -■,)<;- • (33) 

where we used the second line of (1321) . It is straightforward to calculate the derivative of 

d _ 4a g eo(e - |) 2 + a g e ^ + (1 - a g ) 

—n{e — -— . [64: 

de (1 - a g ){\ - e y 

Since this is manifestly positive in the range < eo < 1, we have the inequality h(eo) > 
h(0) = 3, that is, 

m 2 eS > 3H 2 . (35) 

The effective mass of the massive graviton is bounded by Hubble scale from below. Using 
the analysis in Sec. Ill B\ we can see that the anisotropy rapidly decays in a Hubble time. 



IV. DECAY OF THE ANISOTROPY: CASES X f / 

In this section, we construct de Sitter solutions with Xf ^ 0. Then, we check the pertur- 
bative stability of the de Sitter solutions. Finally, we examine if the effective mass of the 
massive graviton can be smaller than Hubble scale. 
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A. de Sitter solutions 



We study de Sitter solutions and give a classification of them. What we should check is 
whether roots of g(e) = are positive and satisfy Hq > 0. 

1. When are roots of g{e) = positive? 

Since the behavior of g(e) is largely determined by the leading term, A/ + C,a g , we discuss 
the following three cases separately. 

1. In the case \f > —£,a g , the coefficient of the leading term in g(e) is positive, which 
indicates 

g(e) —> — oo as e — > — oo , g(e) — > +oo as e — > +oo . 

Combining the above with g(0) — £(1 — a g ) > 0, we see that there always exists 
a negative root. Since g"(0) = — 6£a 5 < 0, the inflection point must exist in the 
positive side of e. Therefore, the number of positive solutions can be characterized 
by the discriminant of g(e) = 0. If the discriminant is zero, a multiple positive root 
exists. On the other hand, if the discriminant is positive, two positive roots exist. The 
discriminant of g(e) — is given by 

D = -27(1 -a g ) 2 (^ + a g y 

+2c[2c 2 + 27a 9 (l - a g )] + a 9 ) + 9a 2 g [c 2 + 12a 9 (l - a g )\ , (36) 

where we defined 

c = \ g /i - 2a g + (1 - a g ) . (37) 
The condition that the discriminant is non negative reads 

A- < A/ < A+ , (38) 

where we defined 

T + a ° = 27(1 - a g )^ 5[2c?2 + 27S(1 " aa)] ± + 9 ° 9(1 " a * )]§ } • (39) 
We can see A_ < — £a 5 and A + > — £a g from (1391) taking into account the inequality 

|2[c 2 + 9a 9 (l-a g )]5| - \c[2c 2 + 27a 9 (l - a g )} \ > . 
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Thus, for A/ = A + , there exists a single multiple positive root of g(e) = 0. Since we are 
considering the range A/ > —£,a g , there exist two positive roots for — £a 9 < \f < A + . 

2. In the case A/ = —£,a g , g(e) becomes the quadratic function of e. Since the coefficient 
of the leading term — 3£a 9 is negative and g(0) = £(1 — a g ) > 0, there exists a single 
positive root. 

3. In the case A/ < — £a 5 , the coefficient of the leading term in g(e) is negative, which 
leads to 



Because of the fact g(0) = £(1 — a a ) > 0, there always exists a positive root. Since 
g"(0) = — 6£a 5 < 0, the inflection point exists on the negative side of e in this case. 
Thus, other possible roots should be negative. Namely, there exists a single positive 
root for A j < — £a 9 . 

We found that two positive roots exist for — £a g < A/ < A + and a single positive root 
exists for A/ < — £a g and Xf = A+. 

Next, we check whether these roots satisfy the condition Hq > 0. 

2. Is Hi > satisfied ? 
Rewriting the first line of ( f!5l) as 



we see that A 9 > — £a fl /4 is at least needed for > 0. Therefore we assume A 9 > — £%/4 
below. Then, we can factorize f l4"0l) as 



g(e) — > +oo as e — > — oo 



(7(e) — 00 as e — )• +00 . 



Hi = A 9 + £a 5 (2-e )(e -l) 




(40) 



Hi = -£a g (eo ~ e P )( e o - e m ) , 



(41) 



where we defined 




(42) 
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Note that e p and e m do not depend on A/. Thus, in order to have Hq > 0, we have to seek 
positive roots of g(e) = in the range 



< e < e 



v ■ 



(43) 



As we discussed in the previous subsection, A/ < A + is needed for the existence of positive 
roots. We first consider the case Xf = A + for which there exists a single positive root. In 
this case, we have to solve g(e*) = g'(e*) = which give rise to the equation 



S e * + 3 ce * - I 1 - a g) = 

The positive root of this equation is given by 

-c+ v/c 2 + 9a g (l 



. 



3a n 



(44) 



(45) 



Thus, we see 



2£(^ + Wc 2 + %(1 



:t + T) + K 



a q ) + yjc 2 + 9a g (l - 



> 



(46) 



9 ' 9 ^ 9 s 

Therefore, the inequality e m < e* < e p must hold. 

As we decrease Xf with fixing A 9 ,a 3 ,£, the discriminant of g(e) = becomes positive. 
Thus, there will be two positive roots until A/ reaches —£,a g . We shall call smaller one inner 
root, e- m , and the other one outer root, e out . We note that e in is always smaller than e* and 
e out is always larger than e* because 

g(0) = f(l - a g ) > 0, p(e„) = ^(A/ - A+) < 0, 5 (e)4+oo as e^+oo. 

We can regard Xf < —£,a g the inner root because the inner root is continuously 

connected to the positive root for A/ < — £a g when A/ crosses — £a g below. 

We shall evaluate the first derivative of €q with respect to Xf since we want to know the 
behavior of the roots when we decrease A/. Differentiating g(Xf, e (A/)) = with respect to 

A/ 

dg{X f ,x(X f )) 



dAy 







x=e 



we obtain 



dfo 
dA, 



dffQz) I 
drr I x=eo 



(47) 



(48) 
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FIG. 1: We plotted g(e) for X g > 2£a g . We set a g = 0.5, A 9 /£ = 1.1. Then A+/£ ~ 1.41. As A/ 
decreases, the outer root increases and the inner root decreases. When Xf reaches X p , the outer root 
crosses e p above and i?o( e out) becomes negative. But the inner root always satisfies HQ(e- m ) > 
since e m is non positive. 



First, we discuss the outer root. Since g(e*) < and g(e) — > +oo as e — > +oo, the outer 
root always satisfies 

> . (49) 



dg(x) 



dx 



X — Gout 



Then, from (HH1), we can see 



de out 



< . (50) 



dA/ 

Therefore, e ou t starts from e* at Xf = A + and e out monotonically increases as A/ decreases. 
We can expect that e out sometime reaches e p . Indeed, e out reaches e p when Xf becomes small 
as 

A p = £(l-a 9 )^^>0, (51) 



where we used the fact Hq = at e p . Therefore, e out exists in the range (e m , e p ) if and only 
if Xf > X p . We mention that A p — > +0 when X g — > +oo since e p — > +oo (see f l42|) ). 

Next, we discuss the inner root. In turn, since g(0) = ^(1 — a g ) > and g(e*) < 0, the 
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FIG. 2: We plotted g(e) for -±£a s < X g < 2£a g . We set a g = 0.5, X g /£ = 0.5. Then A + /£ = 0.5. 
As A/ decreases, the outer root increases and the inner root decreases. When A/ reaches X p , the 
outer root crosses e p above and i?Q(e ou t) becomes negative. When A/ reaches A m , the inner root 
crosses e m below and HQ(e in ) becomes negative. 



inner root always satisfies 



Then from (148|) . we can see 



dg{x) 



dx 
de, 



< . (52) 



dA/ > . (53) 

Therefore, e m starts from e* at Xf = A + and monotonically decreases as A/ decreases. Note 
that e in ( ii ^) 1/3 -> +0 as X f — > — oo. We can expect that e m sometime reaches e m . 
To see this, we need to notice that 

tm = g-./^7i = g-./^- 2 ^ + g (54) 
2 Vf«j 4 2 V 4 

changes the sign at A 9 = 2£a g . Hence, we can consider the following two cases. 

1. In the case X g > 2^a g , e m is non positive. Then e in cannot reach e m when we decrease 
Xf. Therefore, e in always exists in the range (e m , e p ) and satisfies Hl{e- in ) > (see Fig. 

hd. 
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2. In the case — i£a g < X g < 2£a 9 , e m is positive. Then e in can reach e m when we decrease 
Xf. Indeed, e in reaches e m when Xf becomes small as 

A m = £(1 - a 3 ) £m 3 1 , (55) 

where we used the fact Hq = at e m . Therefore, ei n exists in the range (e m , e p ) 
and satisfies Hq > if and only if Xf > X m . We mention that A m — > — oo when 
A g — > 2£a 3 — because e m — >■ +0 (see ( )54j) ). In Fig. EJ we illustrate these features. 

We note that X p > X m when — \^a g < X g < 2£a g . We can see this from the definitions of 
X p and X m as 

8(f^ + if 2 
( 2 ~ jt g ) 

We summarize the results derived in this subsection in Table [H Table [TT] and Fig. [3j 
TABLE I: For X g > 2£a g TABLE II: For -\ta g < X g < 2£a g 





inner 


outer 






inner 


outer 


x + <x f 


X 


X 




A + < Xf 


X 


X 


x f = x + 


o 




X f = A+ 


o 


Xp < Xf < A+ 


o 


o 




X p < Xf < A+ 


o 


o 


A/< A p 


o 


X 




A-m < Xf < X p 


o 


X 








Xf < x m 


X 


X 



In the tables, "o" means there exists a positive root of g(e) = which satisfies Hq > 0, 
i.e., a de Sitter solution exists. And, "x" means there exists no positive root or there exists 
a positive root for g(e) = but Hq < 0, i.e., no de Sitter solution exists. For X g < —^a g , 
there is no root satisfying Hq > 0. Surprisingly, we have an upper bound for Xf and there 
exist de Sitter solutions even for arbitrary large negative Xf in the case X g > 2£a g . We note 
that X g — 2£a 3 can be interpreted as an effective cosmo logical constant if we see the explicit 
constant term in the first line of ( JTBl . It is remarkable that there also exists a de Sitter 
solution for the case effective cosmological constant is zero. 
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FIG. 3: We depicted the region de Sitter solutions exist. We chose a g = 0.5. The multiple solution 
exists on Xf = A+ curve. The outer root exists in the region below A+ and above X p . The inner 
root exists in the region below A + and above X m . X p — > +0 as X g — > +oo and X m — > — oo as 
X g — > 2£a g — as we mentioned in the text. The triple point is given by (X g , A/) =(— |£a 9 , 
^£(1 - a g )) and there Hq = 0. 

B. Stability of de Sitter solutions 

In this subsection, we examine the stability of de Sitter solutions. In Sec. Ill Bl We saw 
that the sign of determines the stability of de Sitter solutions, i.e., solutions are stable 
if m^ e is positive and unstable if is negative. Recalling the formula 

m cfi = £ % e o + (1 - a g )— (3 - 2e ) , 

we can see that m 2 ^ is positive when e < | and m^ s is negative when €q > |. From now 
on, we suppose X g > —j£a g so that Hq > is satisfied. 

We know that g(e) = has positive roots when Xf < X + . We first consider Xf = X + 
case where there exists a multiple positive root e*. Since we supposed X g > —\£,a> g , we can 
evaluate e* as 

3 _ £ = 3(^ + %) > q 

2 "* (^ + a f) + \a g + (1 - a g ) + v/c 2 + 9a g (l - a g ) > 

16 



Then, is positive. Therefore, we find that the de Sitter solution corresponding to the 
multiple root is stable. 

Next, we decrease Xf from A + . The inner root always satisfies e m < e* as we mentioned in 
subsection IIV A 2l We know e* is smaller than |. Therefore the inner root is always stable 
since 6i n < e* < |. On the other hand, the outer root always satisfies e out > e*. Since e* 
is smaller than | and the outer root monotonically increase as A/ decreases, we can expect 
that e out sometime reaches |. Once e out reaches |, the effective mass vanishes. There, A/ is 
given by 



A 



1 



3 



2 27 L V y 4 

and Hubble scale reads 



3(A ff + ^)+e(l-a 9 ) >0, (58) 



HtQ=\ g + ^L>0. (59) 

Note that A p < A| < A + because e* < § < e p (see f|4"2]) and fl57]) ). Therefore, the outer root 
is stable when Xf > A a and unstable when Xf < X$. 

1 2 J 2 



C. Appearance of the Higuchi bound 

In this subsection, we will evaluate the effective mass of the massive graviton correspond- 
ing to the anisotropy. 

From the definition of m^ ff and the first line of ffTBI) . we can deduce the following expression 

2 



m e 2 ff (e ) - 2tf 2 (e ) = 

eo 

3£a 9 



O 9 £o + g Ce — (1 — Clg; 



e* - e )(e - e 2 ) , (60) 



a 9) 



where e* is given in ( 145]) and we defined 

e2 = -c-^+iMr _„ <o 

Since e2 is negative, the sign of m 2 ff — 2Hq depends on that of (e* — eo). Namely, eo = e* 
is equivalent to m 2 ff = 2Hq, e < e„ leads to m 2 ff > 2Hq, and e > e* leads to m 2 ff < 2/Jq. 
When A/- = A + , the multiple root e* obviously satisfies = 2Hq. When A/ < A + , there 
are two positive roots for g(e) = 0. The inner root always satisfies e in < e* as we mentioned 
in subsection IIV A 21 Hence, the inner root always satisfies m 2 ff > 2Hq . On the other 
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hand, the outer root always satisfies e ou t > e*- Therefore, the outer root always satisfies 

m e 2 ff < ml 

Remarkably, the equation m^ s — 2Hq = coincides with the equation determining the 
multiple root e* (see (|44p and (|60p ). That is the reason why the bifurcation point of de 
Sitter solutions is exactly the same as the Higuchi bound. 

Note that the anisotropy decays more rapidly than Hubble time scale 1/Hq for the inner 
root and it decays more slowly than 1 /H or exponentially grows for the outer root if we 
use the analysis of Sec. Ill Bt 

Finally, we shall see that the ratio of the effective mass to Hubble scale monotonically 
varies along the line that the value of eo is constant on X g -Xf plane. We define ( as the ratio 
of the effective mass to Hubble scale, 

t=™k = ^ Q g 6 o + ( 1 - Q g)^K3 - 26 ) 

C HI A 9 + £a g (2-e )(e -l) ' 1 j 

From this expression, it is obvious that <9C/<9A g | eo=const . < for e > § where ( is positive, 
and d(/d\ g \ eo=const , > for e < § where ( is negative. 

We will check how the line that eo is constant can be drawn on X g -Xf plane. When we 
fix the value of eo, g(eo) = gives the relation between X g and A/ as 

X f = \x g - ia g + 3^ + (l-3a^o-e(l-a g ) (g2) 

On X g -Xf plane, each point in the region A/ < A + determines two lines: one for the inner 
root and the other for the outer root. From the fact that 



dA + _ / c+y/c* + 9a g (l-a g ) \ 2 _ 1_ 
d^ _ V 3(1 - a 9 ) J = el 



(63) 



each line is tangential to A/ = A + curve. We also know that A/ = A + is a convex function 
since from 

tt 1 = ; ~ € * = < . (64) 

dA 9 £y/c? + 9a g (l-a g ) 

we can obtain 

d 2 A+ -2de* 

w = t5a; >0 ' (65) 

Therefore, the series of lines cover the whole region satisfying A/ < A+. 
Using these formulae, we depicted Fig. @] and Fig. [5] 
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FIG. 4: We plotted £ =const. curves of the inner root on X g -Xf plane. We set a g = 0.5. The 
inner root is always stable since is positive. In this figure, ( = 2 and C = +oo curves coincide 
with Xf = A+ and A/ = A m curves, respectively. Note that if we start from a point on Aj = A+, ( 
monotonically increases along e; n = const, line. 

V. CONCLUSION 

We investigated the cosmic no-hair conjecture in the presence of spin-2 matter. More 
precisely, we studied the cosmic no-hair conjecture in bimetric gravity using the perturbative 
method. First, we analyzed de Sitter solutions and found that there are two branches of 
de Sitter solutions. We examined the stability of de Sitter solutions and found that there 
always at least one stable solution. Finally, we evaluated the effective mass and found that 
the stable branch of de Sitter solutions satisfies the Higuchi bound. The other branch does 
not satisfy the Higuchi bound. The bifurcation point of two branches exactly coincides with 
the Higuchi bound. Thus, we concluded that there exists a de Sitter solution for which the 
anisotropy decays and the effective mass for these perturbations satisfy the Higuchi bound. 
Since the cosmic no-hair conjecture is already violated at the linear level in known cases, 
our result indicates that the cosmic no-hair conjecture is correct in bimetric gravity even 
though we have not given the nonlinear analysis. 
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1 2 

FIG. 5: We plotted ( =const. curves of the outer root on A 9 -A/ plane. We set a g = 0.5. The outer 
root is stable above A3 where m^ s is positive and unstable below A3 where m^ s is negative. In 
this figure, ( = 2, £ = and ( = —00 curves coincide with A/ = A+, Xj = A3 and Aj = A p curves, 
respectively. We see that if we start from a point on A/ = A+, C monotonically decreases along 
fout = const, line in the stable region and monotonically increases in the unstable region. 

As a future work, it would be interesting to explore the meaning behind the curious fact 
that the bifurcation point of two branches of de Sitter solutions coincides with the Higuchi 
bound. Moreover, since we have found that at least one branch of de Sitter solutions is 
stable in our analysis, we can consider inflation in bimetric gravity without pathologies. It 
would be important to clarify what kind of signatures peculiar to bimetric gravity appear 
for example in the cosmic microwave background radiation. 
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Appendix A: Derivation of Basic equations 

In this appendix, we derive basic equations. 

1. Ansatz and Lagrangian 



We start with the anisotropic metric ansatz for g^ u and f t 



ds i = - N ^ t )dt 2 + e 2a ^[e~ 4 ^dx 2 + e 2 ^\dy 2 + dz 2 )} 



and 



ds 12 = -M 2 (t)dt 2 + e 2 ^[e-^dx 2 + e 2X ^\dy 2 + dz 2 )} , 
respectively. From these metrics, scalar curvatures are calculated as 



Moreover, g 1 f is given by 
/ (M/N) 2 



\ 



3 2^-2«-4A+4cr 



3 2/3-2o+2A-2ct 



3 2/3-2a+2A-2<7 



/7 2 



A 2 



(Al) 



B 2 ) 



where we have defined variables as 



7 = M/N , e = e?- a , rj = e x - a , 



A = e V ~ 2 = e ?- a - 2X+2 ° , B 



erj = e 



/3-a+A-cr 
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Thus, we obtain 



7 



I - A 



1 - B 



V 

Then, we can calculate the interaction term as 

F2 = \{[Lf-[L 2 }) 
1 



1-Bj 



[(A-A-2B- 7 ) 2 - (4 - 2A - 4B + A 2 + 2B 2 - 2 7 + 7 2 )] 
[6 - - 65 + B(2A + B)+ 7 (-3 + A + 25)] . 



(A2) 



(A3) 



Therefore, the Lagrangian reads 

3 , 



C 



N 



a 2 + d 2 ^ 



A^A, 



+ Mie 



^t(-/3 2 + A 5 



MAi 



+m 2 M 2 Ne 3a [Q-3A-QB + B(2A + B) + 7 (-3 + A + 25 



(A4) 



2. Equations of motion and constraints 



We normalize parameters and time with M e as follows: 



Ml 



m 
M 2 



A, 



A, 



3M e 2 



3M e 2 



•/M e . 



Note that < a g < 1 from the definition of M e . 



From the Lagrangian, we obtain the equations of motion 
/o/\' a' 2 1 

\NJ + 3 ^V + Q^ 9 W(3A + 65 - 2B(2A + B)) - M(9 - 24 - 45)] = , 



(A5) 



f 3 ^7 " " a 9)\[ N ^ A + 65 " 25 ( 2 ^ + 5 )) " M ( 9 " - 45)] = , (A6) 







< 




A' N 






)' + 3 


M, 





qV 1 



A 3 
0'X' 1 



£a 5 (A - 5)[A(3 - B) - M] = 



£(1 - a 9 )-(A - B)[N(3 - B) - M] = 



M 3^~ ^e 3 



and the constraints 



\ + ^a 9 [-6 + 3A + 6B - B{2A + B)\ 



(A7) 
(A8) 

(A9) 
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It is easy to find the consistency relation 

M /3'(3A + QB- 2B(2A + B)) - \'{2A - 2B)(3 - B) 



N a'(9 -2A- 45) - a'(2A - 2B) 

From the linear combination of ( 1A7|) and ( 1A8|) . we can also obtain a conserved quantity 



(All) 



1 e 3a a' 1 e 3/3 A' 

E = + . (A12) 

a g N 1 — a g M K J 
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